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Abstract. This paper shows how to determine all the unidimensional two-state 
cellular automaton rules of a given number of inputs which conserve the number of 
active sites. These rules have to satisfy a necessary and sufficient condition. If the 
active sites are viewed as cells occupied by identical particles, these cellular automaton 
rules represent evolution operators of systems of identical interacting particles whose 
total number is conserved. Some of these rules, which allow motion in both directions, 
mimic ensembles of one-dimensional pseudo-random walkers. Numerical evidence 
indicates that the corresponding stochastic processes might be non-Gaussian. 
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1. Introduction 

Systems which consist of a large number of simple identical elements evolving in time 
according to simple rules often exhibit a complex behavior as a result of the cooperative 
effect of their components. Cellular automata are models of such systems. They may 
be defined as follows: Let s:ZxNi— >{0,l}bea function that satisfies the equation 

s(i, t + 1) = f{s{i - ri, t),s{i-ri + l,t),..., s{i + r^, t)), (1) 

for all i & Z and all t G N, where Z is the set of all integers and N the set of 
nonnegative integers. Such a discrete dynamical system is a two-state one- dimensional 
cellular automaton (CA). The mapping / : {0, lyn+'^r+i _^ {0, 1} is the rule, and the 
positive integers r; and rr are, respectively, the left and right radius of the rule. / will 
also be called an n-input rule where n — ri + rr + 1- The function Sf : i ^ s{i,t) is the 
state of the CA at time t. S — {0, 1}^ is the state space. An element of the state space 
is also called a configuration. Since the state St+i at time i -|- 1 is entirely determined 
by the state St at time t and the rule /, there exists a unique mapping Ff : S S such 
that St+i — Ff{St). Ff, which is the evolution operator, is also referred to as the global 
CA rule. 

CA have been widely used to model complex systems in which the local character 
of the rule plays an essential role (Wolfram 1983, Farmer et al 1984, Manneville et 
al 1989, Gutowitz 1990, Boccara et al 1993). In the past few years, CAs have been 
successfully used to model highway traffic. One of the simplest model is defined on a 
one-dimensional lattice of L sites with periodic boundary conditions. Each site is either 
occupied by a vehicle, or empty. The velocity of each vehicle is an integer between 
and Vmax- If x{i,t) denotes the position of car i at time t, the position of the next car 
ahead at the same time is x{i + l,t). With this notation, the system evolves according 
to a synchronous rule given by 

x{i,t + l) ^x{i,t) + v{i,t + l), (2) 

where 

v{i, t + 1) — min {x{i + l,t) — x{i, t) — 1, x{i, t) — x{i, t — 1) + a, Vmax) (3) 

is the velocity of car i at time t + 1. x{i + l,t) — x{i, t) — 1 is the gap (number of empty 
sites) between cars i and i -|- 1 at time t, x{i,t) — x{i,t — 1) is the velocity v{i,t) of car 
i at time t, and a is the acceleration, a — 1 corresponds to the deterministic model of 
Nagel and Schreckenberg (1992) while the case a — Vmax has been considered by Pukui 
and Ishibashi (1995). In this last case, the evolution rule can be written 

x{i, t + 1) ^ x{i, t) + min {x{i + l,t) — x{i, t) — 1, ■Umax)- (4) 

This is a cellular automaton rule with, at least, its left radius equal to Vmax and its right 
one equal to t'max — 1- The case a < t'max is a second order rule, that is, the state at 
time t + 1 depends upon the states at times t and t — 1. For i^max = 1, these two rules 
coincide with elementary CA rule 184 (rule code numbers as in Wolfram 1994) 
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Since, for these highway traffic models on a ring (we shall always consider cyclic 
boundary conditions), the number of cars is conserved, it might be interesting to address 
the more general question: Is it possible to determine all one-dimensional two-state CA 
rules which conserve the number of active sites? We cannot expect that all these rules 
will mimic realistic highway traffic. It is preferable to view them as describing the 
evolution of systems which consist of a fixed number of interacting particles. 



2. General considerations 



If the sites are either all inactive or all active, they should remain so during the evolution. 
Therefore, for any number of inputs n, the local rule should satisfy the conditions 

/(0,0,0,...,0) = (5) 

V 

n 

/(1,1,1,...,1) = 1. (6) 

n 

If the rule (|l|) changes the site value s{i,t), we may say that it either "created" a 
particle, if s{i,t + 1) = and s{i,t + 1) = 1, or "annihilated" a particle in the opposite 
case. Since, the argument s(i, t) of function / takes the values and 1 an equal number 
of times, conservation of particles number implies that the number of creations and 
annihilations should be equal. In other words, the number of preimages of and 1 by 
/ should be the same. 

Consider rules fi and /2, whose radii are, respectively, r^, and r^i, and ri2 and rr2- 
The rule /i o /2 which consists, at each time step, in the successive application of /i 
and /2, conserves the number of particles if /i and /2 do. Its radii are r/ = rn + ri2 
and Tr = Tri + rr2- For instance, the 4-input rule whose binary code number is 
1011100010111000 (n = 1, ri = 2) conserves the number of particles since it is the 
composition of the left shift (binary code number 1010, r; = 0, r,, = 1 ) and rule 
184 (binary code number 10111000, r; = = 1) which both conserve the number of 
particles. 

If, as for highway traffic, we wish to follow particles motion, it might be useful 
to define a representation of rule / which exhibits this motion. Such a "motion 
representation" may be defined as follows. List all the neighbourhoods of a given 
particle represented by 1. Then, for each neighbourhood, indicate the displacement 
of this particle by an integer v, where v is positive if the particle moves to the right 
and negative if it moves to the left. For instance, the motion representation of Rule 184 
would be 

10 1, 11 0. (7) 

Since, for this particular rule, the particle can only move to the right, we only need to 
indicate the relevant neighbourhood of the particle. This representation can be made 
more visual if we draw an arrow joining the initial and final positions of the particle, 
i.e., for Rule 184 

10 , 11 . (8) 
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Note that, in this case, it is not necessary to specify the moving particle by a bold digit. 

This last notation is very compact. For instance, the 4-input rule which results 
from the composition of Rule 184 and the left shift, is represented by 

10 , -ll . 

where • represents either or 1. The motion representation has another advantage. 
When we are interested by the motion of the particles, the knowledge of the rule table, 
which gives the images of the various n-inputs, is not sufficient. We have to specify the 
values of the right and left radii since modifying ri and at constant n is equivalent to 
adding a constant velocity to all the particles. 

Rules obtained by reflection or conjugation of a rule conserving the number of active 
sites have the same property. Reflection exchanges the values of ri and and changes 
the sign of the velocity. Conjugation exchanges the roles of Os and Is, that is, if a rule 
describes a specific motion of particles (represented by Is) then its conjugate describes 
the same rule, but for the motion of holes (represented by Os). If R and C denote, 
respectively, these two operators, two n-input rules /i and /2 are said to be equivalent 
if there exists an element g of the four-group generated by R and C which transforms 
/i into /a. 



3. Rules determination 



One method to determine all the n-input rules / conserving the number of active sites 
is to find a system of equations whose solutions are all the functions 

/: {0,1}" ^{0,1} (9) 

which, for all L > n, satisfy the conditions 

/(Xi, X2, . . . , Xn) + f{x2, Xs,..., Xn+l) H h /(xl, Xi, . . . , Xn-l) 

= Xi+X2-\ hXL, (10) 

for all L-ring configurations (cyclic permutations). Such a system will be called an L- 
system of equations. Conditions (|l^) are clearly necessary, but does it exist a minimum 
value Lmin of L such that they are also sufficient? 

We shall prove that Lmin exists, and is equal to 2n — 2. That is, the necessary 
and sufficient condition for a rule f to conserve the number of active sites is to satisfy 
Relations j ^Tdj ) for L = 2n — 2. 

Before giving a formal proof of this result, we shall present a simple, but not 
rigorous, argument. Given the states s(l, t), s(2, t), . . . , s{n, t) of sites 1, 2, . . . , n at time 
t, the state s{rri+i,t + 1) of site r; + 1 at time t + 1 is determined {n = ri + + 1). 
To determine the states of sites 1 and n at time t + 1, we also need to know the states 
at time t of the r; sites on the left of site 1 and the sites on the right of site n. To 
obtain the minimum number of sufficient conditions satisfied by (|^), we shall require 
that the minimum number of sites we have to add to the original n sites should be such 
that their state values at time t + 1 should depend on, at least, one of the site values 
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s(l,t), s{2,t), . . . ,s{n,t). This condition implies that we should consider an Lmin-ring 
in which the sites 1 — r/ and n + r.^ coincide. Therefore, L^[^ = ri + n + — 1, that is, 

Lm'm 277. 2. 

To prove the above result in a more rigorous way, we shall show that, if L > 2n — 2, 
any equation of an L-system is a linear combination of three equations belonging, 
respectively, to (L — 1)-, {2n — 3)-, and {2n — 2)-systems. More precisely, for all L- 
ring configurations {xi, X2, ■ ■ ■ ,xl}, Equation ([ToD can be written 

(^f{Xi,X2, ...,Xn)+ f{x2, Xs, . . . , Xn+l) H h f{xL-l, Xi, . . . , Xn-l)) 

- if i^l 

+ fix2, X3, . . . , Xi-n+s) + ■ ■ ■ + f{Xn-2, ^L-n+l, • • • 5 ^L-l) 

+ f{xL-n+l, XL-n+2, • • • , XL-1, Xi) + . . . + f{xL-l, Xi, . . . , Xn-2, XL-^n+l)) 

, X2, . . . , Xn-2, ^L-n+li XL-n+2) + f{x2, X^, ■ ■ ■ , X^^n+s) + " " " 
+ f{Xn~2, Xi^n+l, • • • 5 Xl^i) + /(x^-n+l, • • • 5 aJ^-i, X^,) + . . . 
+ f{xL-l, XL, Xi, . . . , X„_2) + /(xl, Xi, . . . , X„_2, Xl_„,+i)) 
= (Xi + X2 H h Xl-i) - (Xi H h X„_2 + XL-n+1 H h Xl-i) 

+ (xi H h x„_2 + xl_„+i H \-xl). (11) 

To verify this result, we have to assume that xl_i = xl, which is always the case 
for any cycle, except, when L is even, for the cycle 1010 ... 10. Verifying ( ]TT| ) is then a 
bit tedious but straightforward. By induction, relation ( pA]) shows that any equation of 
an L-system is a linear combination of equations belonging to {2n — 3)-, and {2n — 2)- 
systems. 

The equation corresponding to the cyclic configuration 1010 ... 10 reads 

/(lOlO ... 10) + /(OlOl . . . 01) + ■ ■ ■ + /(OlOl . . . 01) = ^ (12) 



2 



if n is even, and 



/(lOlO ... 01) + /(OlOl . . . 10) + ■ ■ ■ + /(OlOl . . . 10) = ^ (13) 



2 



if n is odd. That is 



/(1010...10) + /(0101...01) = 1 (14) 
if n is even, and 

/(1010...01) + /(0101...10) = 1 (15) 

if n is odd. One of the images by / of the two alternating n-sequences of Os and Is is 
equal to 1, and the other one to 0. 
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4. Examples 

One- and two-ingut rules conserving the number of active sites are trivial. The identity, 
represented by 1 , is the only one-input rule, and the left and right shifts, represented 
respectively •! and !• , are the only two-input rules. Note that the rule represented 
by !• or •! is the identity viewed as a two-input rule, but in agreement with our 
convention to only represent the relevant neighbourhood, we shall always represent it as 
a one-input rule. This is a general feature. When we solve the system of equations (0) 
for n = 3 and L^^^^ = 4, we will re-obtain the identity, and the left- and right-shifts as 

3- input rules. 

4.I. 3-input rules 

The only nontrivial 3-input rules conserving the number of active sites are Rules 184 
and 226, represented respectively by 

10 , 11 , and 01 , 11 . (16) 

Rule 226, which can be obtained either by reflection or conjugation of Rule 184, models 
exactly the same deterministic highway traffic rule. The only difference, clearly shown 
by the motion representation, is that cars move to the right instead of moving to the 
left. 

4- 2. 4-'i''>T'Put rules 

The system of equations (p!0D for n = 4 and L^in = 6 has 22 solutions. Among these, 
we re-obtain the identity, the left- and right-shifts. Rules 184 and 226 and some simple 
combinations of these rules viewed as 4-input rules. The new rules are: 

• Rules 43944, 65026, 59946, 49024. The motion representation of Rule 43944 
{vi = 2, rr = 1) is 

^ ^ Q 

100 , 101 , 11 . 

This rule coincide with the highway traffic rule (^) for v^ax = 2, and cars moving to 
the right. Rule 65026, obtained by reflection of 43944, describes the same highway 
traffic rule but for cars moving in the opposite direction. 

The motion representation of Rule 59946, which is the conjugate of Rule 43944, is 
Oil , 101 , 111 . 

It describes a highway traffic rule in which drivers, anticipating the motion of the 
car ahead, may move to an occupied site with t>max = 1- More general rules of this 
type have been studied by Fuks and Boccara (1997). Rule 49024 is obtained by 
reflection of Rule 59946. 

• Rules 58336, 52930, 63544, 48268. The motion representation of Rule 58336 

(n = l,rr = 2) is 

^ o o 

100 , 101 , 11 . 
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It describes a highway traffic rule of overcautious drivers who move to the right 
with a velocity equal to 1 if, and only if, they have two empty sites ahead of them. 
By reflection we obtain Rule 52930 describing the same highway traffic rule but for 
cars moving in the opposite direction. 

The motion representation of Rule 63544, conjugate of Rule 58336, is 

^ o Q 

Oil , 010 , 11 . 

A particle moves to the left if, and only if, the neighbouring left site is empty, and 
the neighbouring right site is occupied. If the neighbouring left site is occupied the 
particle does not move. As a highway traffic rule, it describes drivers who do not 
like to be followed, and move to an empty site only when there is a car just behind 
them. By reflection we obtain Rule 48268. 

• Rules 56528, 57580, 62660, 51448. The motion representation of Rule 56528 
(n ^'i-,rr^ 2) is 

100 , 101 , 11 . 

A particle moves to the right if, and only if, its flrst right site is empty and its 
second right site is occupied. As a highway traffic rule it describes drivers who 
move to an empty site if, as a result, they can be just behind another car. Rule 
57580 is obtained by reflection. 

The motion representation of Rule 62660, conjugate of Rule 56528, is 
010 , Oil . 

The particle moves to an empty site on its left if, and only if, there is an empty site 
on its right. Rule 51448 is obtained by reflection. 

• Rules 60200, 48770. These rules are self-conjugate. The motion representation of 

Rule 60200 (r, = l,r^ = 2) is 

r> o <^ o 
100 , 101 , Oil , 111 . 

A particle moves to the right if, and only if, it has two neighbouring empty sites 
on that side. If only the flrst neighbouring site is empty, it does not move to avoid 
occupying a site close to another particle. If its first right neighbouring site is 
occupied, then the particle moves to the left when that site is empty. The effective 
interaction between these particles is repulsive. Rule 48770, obtained by reflection, 
describes a similar evolution rule. 

These last two rules have interesting properties. Starting from a random initial 
configuration, after a maximum number of time steps equal to A^/2, where is the 
number of sites, the system evolves on its limit set. This limit set has a rather simple 
structure. If the density of particles p = |, it consists of 3 types of periodic sequences, 
namely: 

. . . 101010101010 . . . , of period 2 
...100100100100..., of period 3 
...110110110110..., of period 3. 
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The probabilities of the various 3-blocks have been determined numerically. We have 
found 

P(OOO) = P(lll) = 0, 

P(OOl) = P(llO) = P(IOO) = P(Oll) = 0.145 ± 0, 001 
P(OIO) = P(lOl) = 0.210 ± 0.001. 

Regarded as a formal language (Wolfram 1984, Denning et al 1978, Hopfcroft and 
UUmam 1986), such a hmit set is regular. Words in a regular language, on the alphabet 
{0, 1}, are generated by walks through a finite directed graph whose arcs are labeled 
with or 1. Given a finite graph, it is always possible to find an equivalent deterministic 
finite graph, that is, a graph in which no more than one arc of a given label leaves each 
vertex. For Rules 60200 and 48770, the corresponding deterministic graph is represented 
in Figure 1. For p = |, when the CA evolves on its limit set, each particle performs a 




Figure 1. Regular language graph for Rules 60200 and 48770. 

pseudo-random walk. The CAs rules being deterministic, the randomness comes from 
the randomness of the initial configuration. Numerical simulations shows that any 
particle has a probability p = 0.29 to move either to the left or to the right, and a 
probability q = 1 — 2p = 0.42 not to move. Actually this pseudo-random motion is 
periodic in time, the period being equal to A^/2. In the limit set, for a given random 
initial configuration, all particles performs the same pseudo-random walk, with a phase 
difference depending on the distance separating them. More precisely, if denotes 
the position of particle n at time t, for Rule 60200, we have 

X,(t)=X„+i(t-l)-2, 

which implies 

This last result shows that the position of a specific particle at time t is determined by 
the position of another specific particle in the initial configuration. 

To characterize the nature of the randomness of the motion of a particle, we have 
determined the Hurst exponent (Hurst 1951, Hurst et al 1965, Feder 1988) of the time 
series generated by the displacement of a given particle. Given a time series s{t), the 
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Figure 2. Fundamental diagram for Rule 60200. Small circles represent numerical 
results. The piecewise linear line has been obtained using local structure approximation 
(see below). 



Hurst exponent H characterizes the asymptotic behaviour of the standard deviation of 
s{t) as a function of time. A Brownian motion (symmetric random walk) has a Hurst 
exponent H — ^. For a particle moving according to the 4-input rules 60200 and 48770, 
we have found H — 0.63 ± 0.02. Since the pseudo-random motion is periodic in time 
with a period equal to half the lattice size, this numerical result is arguable. It could be 
interesting to perform a detailed study of the correlations, but, even correlated random 
walks may have a Gaussian behavior when the number of time steps goes to infinity 
(Weiss 1994). However, in this case, there exists a crossover between a non-Gaussian 
and a Gaussian behavior. This fact implies that, for a large value of the number of 
time steps t, the exponent of the standard deviation of the walk could, numerically, be 
different from 1/2. 

When p ^ ^, the limit set consists of the previous periodic sequences and either 
sequences of Os if p < | or sequences of Is if p > |. Prom the motion representation of 
Rule 60200, it follows that the average velocity (v) of the particles is P(IOO) — P(Oll). 
The conjugacy operator changes (v) in —(f) and p in 1 — p. Therefore, the so-called 
"fundamental diagram" of road traffic theory, that is, the graph of the ffow p{v) as a 
function of the density p, has a center of symmetry, namely, the point {p, p{v)) = (|,0) 
(Figure 2). Rule 48770 has identical properties. 
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4-3. 5-input rules 

The number of rules conserving the number of active sites grows very fast with the 
number of inputs. There exist 428 5-input rules conserving the number of active sites. 
Few of them are not new either because they actually depend upon a smaller number of 
inputs or because they are simple composition of rules already obtained. In this section 
we shall just describe the self-conjugate rules.| 

There exists 20 self-conjugate rules. Some, such as the identity and the shifts (left 
and right, simple and double), are trivial. We also re-obtain the two self-conjugate 
4- input rules. Each of them twice depending on which side, left or right, the extra 
input is added. Finally, we are left with 11 new self-conjugate rules. For each rule we 
shall always choose the values of ri and such that the condition {v) = for p = 1 is 
satisfied. This can always be done. 

Few of these rules are still not very interesting. After few time steps, for all p g]0, 1[, 
5 rules emulate the identity, which means that no particles are moving. These rules are: 
Rule 3464560268 (r^ = 1, = 3), whose motion representation is 

O O /~\ O /\ o o 

0001 , 0111 , 1001 , 0101 , 0110 , 111 , 1101 , 
Rule 3771264248 {vi = 1, = 3), whose motion representation is 

O /\ Q Q Q Q 

0011 , 0010 , 0011 , 0101 , 1101 , 111 , 1011 , 
and Rule 3824738360 (r^ = = 2), whose motion representation is 

/~\ O O Q Q n, Q 

0011 , 0010 , 1101 , 0111 , 101 , 1100 , 1111 . 

Rules 4249668928 and 415766320 obtained by reflection of the first two rules have the 
same property. Rule 3824738360 is invariant under reflection. 

Rule 3167653058 (n = 3, = 1), whose motion representation is 

^ o o o o ^ 

0001 , Oil , 101 , 1111 , 1001 , 1110 , 

is rather peculiar. As shown in Figure 3, this rule emulates the identity only for p G [|, |]. 
Rule 4270014080, obtained by reflection has identical properties. Note that the flow 
diagram (Figure 3) is piecewise linear. 

The 4 remaining rules are similar to the 4-input rules 60200 and 48770 in the sense 
that they have similar flow diagrams and that, for p = 1, they mimic pseudo-random 
walkers. These rules are: 

• Rule 3221127170 (r/ = 2, = 2), whose motion representation is 

0001 , 1111 , 1001 , 101 , lib , ino , 

and Rule 3937086120 obtained by reflection. 

• Rule 3707031748 (r/ = 2, = 2), whose motion representation is 

,r\ Q Q Q Q /\ Q 

0010 , 0011 , 1100 , 0111 , 101 , 1101 , 1111 , 

and Rule 416291200 obtained by reflection. 
I Codes of all other rules can be obtained from the authors through e-mail. 
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Figure 3. Fundamental diagram for Rule 3167653058. Small circles represent 
numerical results. 
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Figure 4. Fundamental diagram for Rule 3221127170. Small circles represent 
numerical results. 
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Figure 5. Fundamental diagram for Rule 3707031748. Small circles represent 
numerical results. 

The fundamental diagrams of Rules 3221127170 and 3707031748 are represented, 
respectively, in Figures 4 and 5. Here again we verified that the corresponding stochastic 
processes are not Gaussian. We have found that their Hurst exponents are equal for all 
of them to 0.57 ± 0.02. We have no explanation why Rules 3221127170 and 3707031748 
should have the same exponent. 

5. Approximate methods 

The mean-field approximation, which neglects correlations in space and time, yields, 
for these systems, an exact but trivial result. Let p{t) denotes the particles density at 
time t. To find the expression of p{t -|- 1) as a function of p(t), we have to find all the 
preimages of 1 by the n-input rule /. According to (|^) all these preimages contain at 
least one 1. Moreover, among all the preimages containing exactly k + 1 times the digit 1 
{0 < k < n — 1), according to the conditions (0) for L = n, only ('^^^j have a preimage 
equal to 1. Therefore, 

p{t + i) = Pit) ( E l ^) {p{t))\i - p(t))"-^-^) 

= p{t){p{t) + {I - p{t))r' 
=pit), 

which expresses that density is conserved. 

There exists a variety of other approximate methods which, taking into account 
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short-range correlations, improve the prediction of the mean-field approximation. 
Instead of expressing the evolution of the CA in terms of 1-block probabilities, they 
express it in terms of n-block probabilities (Gutowitz et al, 1987). For example, in the 
case of a 4-input rules, the evolution of the 2-block probability distribution is given by 

P{aia2) = w(aia2 \ feoi'i^'2i'3^'4)-P(&o&i&2&3&4), 

where P{aia2) is the probability of block aia2, and 

w{aia2 I bobib2b3b4} = w{ai \ bQbib2b^)w{a2 \ bib2b^bi) 

is the conditional probability that the 4-input rule maps the 5-block 6o^'i^'2&3&4 into 
the 2-block aia2. This equation is exact. The approximation consists in replacing the 
5-block probability P{bQb\b2b^b:i) in terms of 2-block probabilities. That is, 

p., , , , , ^ P{bob,)P{h,h2)P{h2h,)P{h,h,) 

P(bobib2bsbA) = -, w w ^- 

(P(6iO) + P(6il)) (P(620) + P(62l)) (P(&30) + P(&3l)) 

We applied this method up to approximation of order 3 (mean-field being order 1) to 
4-input Rule 60200. The results are not exact, but for the flow diagram the agreement 
with our numerical results is extremely good (Figure 2). 



6. Conclusion 



We have established necessary and sufficient conditions to be satisfied by any one- 
dimensional cellular automaton rule conserving the number of active sites. This result 
has been used to determine all the 4- and 5-input one-dimensional cellular automaton 
rules having this property. These rules express the evolution of one-dimensional systems 
of interacting particles whose number is conserved. Simple deterministic highway 
traffic rules belong to that class of rules. These rules are a natural generalization of 
deterministic traffic rules already studied. We have studied in more detail (flow diagram, 
local structure approximation) some of our rules allowing motion of the particles in both 
directions. When the particle density is equal to |, these rules mimic the evolution of 
pseudo-random walkers. Numerical evidence seems to indicate that the motion of these 
walkers might be non-Gaussian. 
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